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ABSTRACT 



We use a long (300 ksec), continuous Suzakii X-ray observation of the active nucleus in NGC 1365 to investigate the structure of the 
circumnuclear BLR clouds through their occultation of the X-ray source. The variations of the absorbing column density and of the 
covering factor indicate that the clouds surrounding the black hole are far from having a spherical geometry (as sometimes assumed), 
instead they have a strongly elongated and cometary shape, with a dense head (n ~ 10" cm"') and an expanding, dissolving tail. We 
infer that the cometary tails must be longer than a few times 10'^ cm and their opening angle must be smaller than a few degrees. 
We suggest that the cometary shape may be a common feature of BLR clouds in general, but which has been difficult to recognize 
observationally so far. The cometary shape may originate from shocks and hydrodynamical instabilities generated by the supersonic 
motion of the BLR clouds into the intracloud medium. As a consequence of the mass loss into their tail, we infer that the BLR clouds 
probably have a lifetime of only a few months, implying that they must be continuously replenished. We also find a large, puzzling 
discrepancy (two orders of magnitude) between the mass of the BLR inferred from the properties of the absorbing clouds and the 
mass of the BLR inferred from photoionization models; we discuss the possible solutions to this discrepancy. 

Key words. Galaxies; Individual: NGC 1365 - Galaxies: Seyfert — Galaxies: nuclei; X-rays: galaxies 



1. Introduction 

The variability of the absorbing gaseous column density Nh 
observed in the X-ray spectra of Active Galactic Nuclei have 
revealed that a significant fraction of the absorbing medium 
must be clumpy dRisaliti et all 120021) . Moreover, Nh variations 

' on time scales as short as a few days or hours have shown 
that a significant fraction of such absorbing clouds must be lo- 
cated very cl o se to the nuclear X-ra y sourc e (Elvis et al., 2002 

', iRisaUtietail l2005al: iPuccetti et aTL l2007h and, more specifi- 

• cally, within the Broad Line Region (BLR). 

■ The nucleus of the galaxy NGC 1365 is amongst the sources 
that were investigated more thoroughly in this respect, thanks to 
its brightness and probably also to a fortunate (intermediate) in- 
clination of the absorbing medium relative to our line of sight. 

. This is a nearby (z=0.0055) type 1.8 Seyfert galaxy that was 
observed several times in the X-rays, and displayed strong col- 
umn density variations on t ime scales as short as a few days (e.g. 
iRisaliti et an.l2005all2007h . Most of the observations had probed 
absorption variations of this source in discrete, non-contiguous 
time intervals. More recently, (nearly) continuous XMM obser- 
vation s of NGC 13 65 were presented in Risaliti et al. (20093) 
and in IRisaliti et alJ (l2009b). In particular, the latter observation 
allowed the continuous monitoring of an absorption event last- 
ing about 40 ksec, which was modelled with a cloud partially 
eclipsing the X-ray source. More specifically, the spectral varia- 
tions were modelled by assuming a constant column density of 
the absorber (Nh ~ 3.5 10^-* cm"^) and a variable covering fac- 
tor The inferred density, velocity and distance from the black 



hole of the variable absorber is strikingly similar to those of 
the BLR clouds. The complex spectrum always requires a fore- 
ground non-variable (or very slowly variable) absorber (Nh ~ a 
few times 10^^ cm^^), totally covering the source, possibly asso- 
ciated with gas at larger distances from the nucleus (possibly in 
the host galaxy). 

NGC 136 5 has also been observed with Suzaku 
dMitsuda et al.L l2007h with a continuous observation last- 
ing more th an 300 ksec. The integrated spectrum was akeady 
presented in IRisaliti et al.l (l2009c), mostly focusing on the high 
energy part obtained with the Hard X-ray Detector (HXD), and 
revealing an unusually strong excess of emission at E > 10 keV. 
Such excess is ascribed to X-ray radiation piercing a partial, 
Compton thick (Nh ~ 4 x 10^** cm"^) absorber The high ener 
excess appears stable over long time scales (Risaliti et al.. 200i 
I2009cl) . The emerging scenario is that three distinct absorbers 
are present: (1) a distant absorber with Nh < 10^^ cm"^, 
probably associated with gas in the host galaxy; (2) an absorber 
made of broad line region (BLR) clouds with Nh ~ lO^-'^^'^cm"^ 
rapidly orbiting around the black hole; (3) an absorber with 
Nh ~ 4 X 10^"* cm"^ partially covering the source, consisting 
either of the outer region of a warped accretion disk, or of a 
large number of small (< lO'^cm) and dense (n > lO'^cm"-') 
clouds located in the vicinity of the accretion disk. 

In this paper we focus on the absorption variations in the 
spectral range 2-10 keV, ascribed to BLR clouds passing along 
the line of sight, through a detailed temporal analysis of the spec- 
tra provided by the Suzaku low energy (XIS) data of NGC 1365. 
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Fig.l. Merged XISO and XIS3 spectra (1-10 keV) and HXD 
spectrum (15-60 keV) of NGC 1365, integrated over the whole 
duration of the observation. No spectral fit is shown since, as 
discussed in the text, the absorption components are highly vari- 
able, hence a fit (with constant components) to the total, time- 
averaged spectrum is meaningless. 



We identify two eclipsing events that we can model in detail in 
terms of variations of both covering factor and column density. 
The time variations of these quantities are highly asymmetric 
and strongly suggest that the absorbing clouds have a cometary 
shape, i.e. a well defined head and a more diffuse elongated, con- 
ical, tail. 



2. Observations and data reduction 

The obse rvations and the basic data reduction were already re- 
ported in iRisaliti et al.l (l2009cl) . In this section we only briefly 
summarize the observation s and data reduction of the X-ray 
Imaging Spectrometer (XIS i Kovama et al] , l2007h data. 

The observations were performed on 2007 January 21-25, 
for an elapsed time of about 320 ksec and a net exposure of 
160 ksec. The spectra and calibrations were obtained follow- 
ing the standard procedure described in the Suzaku reduction 
guide. The XIS spectra and light curves were extracted in a re- 
gion with a radius of 2.9 arcmin (i.e. matching the optimized 
aperture for the arf matrix) around the nucleus of NGC 1365. 
The background was extracted in a region far from the source, 
with the same radius adopted for the source extraction. The cal- 
ibration files were produced by using the FTOOLS 6.6 pack- 
age. The reduction of the Hard X-ray Detection (HXD) data, 
only m arginally used in this paper, is discussed in IRisaliti et al.l 
(I2009d) . 



3. Temporal and spectral analysis 

In this paper the spectra from the three CCDs (XISO, XIS 1 and 
XIS3) are always treated separately, especially in the statistical 
analysis. However, for sake of clarity, some of the figures will 
show the merged spectra or light curves, for illustrative purposes 
only. 

In the 0.5-10 keV range during the whole observation the 
three XIS collected a total of about 2 10^ photons. Fig.[T]shows 
the XISO and XIS3 merged spectrum integrated over the whole 
observation along with the HXD spectrum. No spectral fit is 



shown for the total, time-averaged spectrum, since the strong 
spectral variability (especially of the absorption components) 
over the entire observation, as discussed in the following, makes 
any fit with constant component s of the integrated s pectrum 
meaningless (as also discussed in IRisaliti et al.Ll2009bl) . 

Fig. 12^ shows the combined light curve in the 2-5 keV band 
rebinned to 2.880 ksec (about half of the Suzaku orbital period). 
The light curve shows strong variations even on short timescales 
of a few ksec, including two sharp drops (marked A and B in 
Fig. |2]i that are strongly reminiscent of the eclipse events seen 
in Chandra and XMM. The egress from these "eclipses" is less 
clear cut however, suggesting more complex structure in the ab- 
sorber. 

To investigate the nature of these variations we have ex- 
tracted the X-ray spectra in 12 time intervals, selected either to 
match the main variations of the light curve or by simply divid- 
ing long time intervals into equally spaced sub-intervals. Within 
each time interval the light curve generally shows additional 
structures, but we decided not to further split the time intervals 
for the extraction of the spectra, since the resulting poorer statis- 
tics would make the errors on the fitting parameters too large, or 
would introduce strong degeneracy between some of them. On 
average, each time interval contains about 1.5 10^ photons (e.g. 
Fig.EJ. 

The resulting 36 spectra (for each time interval there are 3 
spectra, XISO, XISl and XIS3) are simultaneously fitted with 
a model consisting of multiple compon ents, quite similar to 
previous modelling of the same source ( Risaliti et al.L l2009bt 
Guainazzi et al., 2009), with most components and parameters 
held constant for all time intervals (and hence constrained by the 
full statistics of 2 10^ photons) and a few parameters left free to 
vary for the different time intervals. In the following we summa- 
rize the various models components. 

1. Soft thermal compo nents. As reviewed and discussed in 
iGuainazzi et al.l (I2009.) the soft part of the spectrum (0.5- 
2 keV) is ascribed to a complex combination of hot plas- 
mas present in the circumnuclear star forming regions 
(Wang et al., 2009) and in the Narrow Line Region cones, 
as well as the contribution from point s ources (X-ray bina - 
ries and Ultra Luminous X-ray sources, 'Soria et aU 20091) . 
Following the detailed analysis in Guainazzi et al] (l2009t) we 
fit the soft spectrum with three components: two absorbed 
thermal plasmas ("mekal" in XSpec) and a black-body. The 
resulting best fit pa rameters are generally in agreement with 
those obtained bv IGuainazzi et al] (l2009l) (see below), keep- 
ing in mind that the aperture adopted here is different, hence 
we do not expect an exact match. However, the details of 
these soft components are unimportant within the context of 
this paper, since they are constant in time and do not affect 
the variability observed in the 2-10 keV spectral region. 

2. Powerlaw and two cold absorbers. The bulk of the 2-10 keV 
emission is due to a primary powerlaw component, ascribed 
to the corona surrounding the black hole accretion disk, 
with a double absorber: one absorber (Nh,i) totally cov- 
ering the X-ray source, which we shall ascribe to a fore- 
ground medium far from the X-ray source, and a second ab- 
sorber (Nh.2) partially covering the source, which, following 
IRisaliti et al.l y009b), we shall ascribe to clouds transiting 
close to the X-ray source (probably BLR clouds). By using 
wider time intervals (~50 ksec) we found that the spectral 
index does n ot change within the errors (as also found by 
IRisaliti et al.L l2009b). Therefore, we assumed that the spec- 
tral index is the same for all 12 time intervals. The powerlaw 
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normalization is left free to vary among the time intervals. 
The first, total, absorber is kept constant for all time inter- 
vals. The absorber with partial covering is left free to vary 
both in terms of column density (Nh.2) and in terms of cov- 
ering factor (CF). 

3. Narrow FeK line. A narrow (unresolved) iron FeKa line at 
~6.4 keV is clearly detected. In larger temporal bins the flux 
of this line is not observed to vary. As discussed by various 
previous works the narrow FeK line is probably produced by 
circumnuclear dense clouds (most likely Broad Line Clouds 
and circumnuclear absorbing medium in general). The pa- 
rameters of this line are therefore kept constant and fitted 
simultaneously for all time intervals^ 

4. Broad FeK line. As discussed in iRisaliti et all ( l2009bl) the 
spectrum of NGC 1365 (when not in a Compton thick state) 
also requires a broad, relativistic iron FeKa line. The param- 
eters of the iron line (fitted with the "laor" function within 
Xspec) are kept constant for all time intervals, except for its 
normalization that is required to be proportional to the pri- 
mary powerlaw in each time interval (i.e. constant equiva- 
lent width). The physical justification of this assumption is 
that, according to standard models, the relativistically broad- 
ened iron line comes from a very small region, of the or- 
der of a few Schwar zschild radii, (Mbh ~ 2 10^ M0, hence 
Rs ~ 6 10" cm, IRisaliti et all 12009b') and therefore it 
should respond to the continuum variations on timescales 
much shorter (~100 sec) than the size of our time intervals. 
However, in a few representative intervals we have tried to 
leave the normalization of the broad line free, and the fit- 
ting results do not change significantly, the main eff'ect being 
mostly that the confidence intervals increase by about 20- 
30%. 

5. FeXXV and FeXXVI a bsorption lines. As first discovered by 
IRisaliti et all (l2005bh the spectrum of NGC 1365 is charac- 
terized by a set of four absorption lines at ~7 keV due to 
Fe XXV and Fe XXVI, both Ka and Kfi. These lines are 
also clearly detected in our Suzaku spectrum (Fig.[T]i. As in 
IRisaliti et al.l(l2005b l). the depth of each absorption line is left 
free to vary independently in each time interval, but the four 
lines are required to have the same velocity shift within each 
interval. We note that given the high ionization state of this 
absorber, its only spectral components that we expect to be 
important are just the Fe absorption lines, i.e. there is no sig- 
nificant continuum photoelectric absorption associated with 
this absorber. 

6. Reflection components. As in IRisaliti et al.l (l2009bl) the fit 
also includes a set of three (constant in time) continuum 
reflection components: i) a power law with the same spec- 
tral index as the primary continuum, which accounts for the 
free electron scattering by an ionized reflecting medium on 
large scale (e.g. hot gas in the Narrow Line Region), hence 
this component is assumed constant over the whole observa- 
tion; ii) a cold reflection (modelled with "pexrav" in Xspec), 
associated with the same medium responsible for the FeKa 
narrow line emission, likely the circumnuclear dense clouds, 
hence it is assumed to be constant, as in the case of the 
narrow FeK; iii) an inner ionized absorber (modelled with 
"pexriv" in Xspec and subject to the same absorption as the 
primary power law component) to account for warm reflec- 
tion due to ionized gas close to the X-ray source, possibly the 
accretion disk itself; this third component is assumed, as for 
the broad FeK line, to respond to the continuum variations 
on timescales much shorter than our time intervals, hence 
its normalization is assumed to scale proportionally to the 



primary powerlaw continuum within each time interval. We 
will see that the reflection components i) and iii) contribute 
very little to the global fit. Finally, we note that component 
iii) may be relativistically blurred, however given its small 
contribution to the spectrum such "blurring" cannot be con- 
strained and it does not affect the total fit. 
7. Compton thick absorbed powerlaw. As mentioned in sect.[T] 
the Suzaku spectrum at E > 15keV (HXD) shows an ex- 
cess of emission with respect t o the extrapolation of the 
above components Risahti et al.l (l2009d) . This excess is best 
interpreted by assuming that the intrinsic X-ray emission 
is actually stronger than observed in the 2-10 keV spec- 
trum, but partially absorbed by gas very close to the X-ray 
source that we fit with Nh,3 = 3.7^[q IO^'* cm"^. The in- 
ferred covering factor strongly depends on the geo metry of 
the Compton thick medium, as discussed in Risali tiet al.l 
(2009c). Compton thick partial covering has also been ob- 
served in other sources (Turner et al....2009.) . This inner, very 
thick, partial absorber is either associated with gas in a 
warped accretion disk or a population of small and dense 
clouds surrounding the accretion disk. The very high-energy 
excess is found to be stable over long time scales (years), 
hence in our analysis of the Suzaku data this component was 
assumed constant for all time intervals. We note that this ad- 
ditional component at high energies does not aff'ect the fitting 
of Nh,2 and CF. 

The global, simultaneous fit gives a reduced xl-^-^55. 
Overall there are 25 free parameters common to all time intervals 
that, therefore, are constrained by the full 2 10^ photon statistics. 
Each time interval has its own set of 8 free parameters, three for 
the continuum (Nh,2, CF and powerlaw normalization) and five 
for the absorption Fe lines (velocity and individual depth). 

Fig.|3]shows the spectra and the resulting fit in two represen- 
tative time intervals (T8 and T12, as marked in Fig.|2]i, clearly 
displaying a strong variation of absorption between the two time 
intervals. 

In the following we will mostly focus on the variations of 
the column density Nh,2 and of the Covering Factor (CF) of the 
partial absorber, whose resulting best-fit values in the various 
time intervals are shown in Fig.|2]3 and c. 

One possible concern is whether these two quantities are de- 
generate, given the reduced statistics available in each time inter- 
val. We found this not to be the case in any of the time intervals 
(except for the last interval, which does not require a second ab- 
sorber). As an example, Fig.|4|shows the confidence levels in the 
Nh,2 versus CF plane, for the time interval T8: although there is 
some correlation between the two parameters around the best 
solution, the 90% confidence contour provides good constraints 
on both quantities. We emphasize that this is the first time that a 
temporally resolved X-ray spectral analysis is able to break the 
degeneracy between the evolution of the column density (Nh) 
and the covering factor (CF) of the X-ray absorber in an AGN. 

In the rest of this section we report the fitting results for the 
other parameters, which, although important for the global fit, 
are less relevant for the topic if this paper. 

Table [1] lists the best fitting values for the main parameters 
in common to all time intervals. The first part reports the soft 
(not variable) thermal components. The second part reports the 
powerlaw slope and the reflection components. Note that the 
continuum slope F obtained here is different than obtained by 
Risaliti et al. (2009c) on the same set of data, i.e. F = 2.34+°°^. 
There are two clear x^ minima for this parameter, but the de- 
tailed analysis in sub-intervals performed here (hence account- 
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Fig. 2. a) Photon count rate in the soft band 2-5 ke V as a func- 
tion of time (in time intervals 2.88 ksec). The vertical, green 
dashed lines indicate the beginning of the two main eclipses 
produced by two clouds passing in front of the X-ray source. 
The vertical errorbars indicate the Icr fluctuations within each 
bin. b) Column density (Nh,2) of the variable absorber in the 
12 time intervals, resulting from their simultaneous spectral fit. 
Errorbars on Nh.2 are at the 90% confidence. The black dotted 
lines indicate a quadratic fit to Nh,2 in the eclipses "A" and "B". 
c) Covering factor (CF) of the variable absorber. Errorbars on CF 
are at the 90% confidence. The dotted lines show the variation 
of covering factor expected from the cometary model shown in 
Fig. |8] and discussed in the text. The segments labelled with T8 
and T12 show the extraction time intervals of the spectra shown 
in Fig. [3] 

ing for the spectral variability) shows that the x^ minimum with 
F = 1.81 fits the data with much higher significance (xl = 1.05, 
versus ;^'J = 1.08, with about a thousand d.o.f). The reflection 
factor R of the reflection components is given relative to the 2- 
10 powerlaw continuum observed in the last time interval, T12 
(Fig. |2]i. The real value of the reflection factor should be given 
relative to the intrinsic powerlaw component, i.e. the continuum 
observed at E>20 keV corrected for the Compton thick partial 
covering Nh,3, which is much higher than the component ob- 
served in the 2-10 keV range and which gives values of R that 
are much lower However, applying the correction to the contin- 
uum transmitted at high energies is very uncertain, since it de- 
pends strongly on the geometry o f the Compton thick a bsorbing 
medium, as discussed in detail in lRisaliti et al.l (l2009cl) . 
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Fig. 3. Merged XISO and XIS3 spectra of NGC 1365 and result- 
ing fitting model in two representative time intervals, T8 (be- 
tween 167 ksec and 185 ksec, red) and T12 (between 261 ksec 
and 313 ksec, blue), as indicated in Fig. |2] Although, for the sake 
of clarity, we plot the combined XISO and XIS3 spectra in both 
figures, in the fitting analysis we used all individual XIS spectra 
(XISO, XISl and XIS3), kept separate. 




Fig. 4. 68%, 90% and 95% confidence levels in the Nh,2 versus 
Covering Factor diagram, for the time interval T8 (Fig. |2]l. The 
cross marks the best fit parameters. 



The third part of Table \T\ reports the constant absorption 
components. For what concerns the Compton thick absorbing 
medium, Nh,3 only accounts for the photoelectric absorption, 
which is enough for our purposes, since we are not interested 
in the high-energy component. The full treatment including the 
Compton scattering is given in lRisaliti et al.l (l2009d) . 

The fourth part of Table [T] provides the parameters of the 
emission lines. The equivalent widths are given relative to the 
last time interval, T12. 

Fig. |5] shows the variation of the normalization of the pow- 
erlaw continuum in the 2-10 keV range in the same intervals 
analyzed in Fig.|2] The normalization changes significantly dur- 
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Table 1. Best fit parameters common to all time intervals 
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Fig. 5. Variation of the normalization of the powerlaw contin- 
uum in the 2-10 keV range over the 12 time intervals analyzed in 
Sect. |3] The green dashed lines indicate the beginning of eclipses 
A and B, as in Fig.|2l 



ing the observations. This is not unexpected since NGC 1365 
is a Narrow Line Syl and, as such, it is expected to be charac- 
terized by significant variability. It is important to note that the 



Fig. 6. Equivalent width of the Fe absorption lines over the 12 
time intervals analyzed in Sect. [3] The green dashed lines indi- 
cate the beginning of eclipses A and B, as in Fig.|2] 



variations in the normalization of the powerlaw are uncorrected 
with the 2-5 keV light curve behavior, further supporting the idea 
that the variations of the light curve also require variation of the 
absorbing column density and covering factor 

Fig. |6] shows the equivalent width of the four Fe absorption 
Unes (FeXXV Ka 6.697 keV, FeXXVI Ka 6.966 keV, FeXXV 
igS 7.880, FeXXVI K/3 8.268 keV) in the 12 time intervals an- 
alyzed in Sect. [3] Although there is some indication of varia- 
tion of the depth of the absorption lines, especially in association 
with eclipse B, the variations are all within the statistical errors. 
Fig. |7] shows the bulk velocity of the absorption lines, which 
also does not show significant changes within eiTors. We note 
that the infeiTed EW are in line with those obtained for the same 
source by Risaliti et al. (2005b). In the latter work variability of 
the lines (both in terms of EW and of velocity) was observed on 
timescales of a few months. However, our data are not sensitive 
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Fig. 7. Velocity of the Fe absorption lines over the 12 time in- 
tervals analyzed in Sect. [3] The green dashed lines indicate the 
beginning of eclipses A and B, as in Fig.|2] 



to small changes in the absorption lines EW and velocity shift 
on much shorter timescales, with lower exposures. 



4. The cometary shape of the eclipsing clouds 

From Fig. |2] and in particular from the light curve and the Nh,2 
variations, we identify two main eclipsing events, labelled "A" 
and "B", due to the passage of two clouds in front of the X-ray 
source, with eclipse "B" being the more prominent and more 
clear event. The decreasing Nh,2 in the time interval prior to 
ecUpse "A" probably shows the final part of the occultation by 
a previous cloud. The eclipses are far from being symmetric, 
implying that clouds are not spherical. A spherical cloud sig- 
nificantly larger than the X-ray source would either produce a 
sudden total obscuration (CF - 1) or a smooth and symmetri- 
cal evolution of the covering factor if grazing the X-ray source; 
while a spherical cloud with size comparable to, or smaller than 
the X-ray source would simply produce a short, symmetric dip 
in the light curve and in the covering factor 

Instead, both eclipses start with an abrupt drop of the 2- 
5 keV flux (marked by the green vertical dashed lines A and 
B in Fig.|2]i occurring within less than a few thousand seconds. 
The spectral analysis of the first time interval of each eclipse 
(Fig.|2j5 and c) indicates that such a sudden drop of X-ray flux is 
due to the transit of a cloud with Nh ~ 10^^ cm"^ suddenly cov- 
ering about 65% of the source. Subsequently the covering factor 
increases, but more slowly, on time scales of about 50 ksec. 

These observational results indicate that the occulting cloud 
must be elongated, with a size perpendicular to the direction of 
motion comparable to the X-ray source, while the size along the 
direction of motion is at least an order of magnitude larger This 
is an immediate consequence of the fact that the beginning of the 
eclipse is very sharp (with the covering factor jumping to 65% 
in about a ksec) while the subsequent evolution of the covering 
factor is much smoother on timescales of several 10 ksec. 

The fact that the covering factor slowly increases, reaching 
unity about 50 ksec after the beginning of the eclipse, indicates 
not only that the clouds are elongated, but that they must have a 
conical, cometary shape, with their tail gradually occulting the 
whole X-ray source. 

The decrease of column density with time within each 
eclipse indicates that the head of the "comet" is much denser 
than its tail, while the latter fades gradually into a low density 
medium. We note that even the time interval before eclipse "A" 
shows a similar trend, i.e. decreasing Nh,2 and increasing cov- 
ering factor, hence in these intervals we may be observing the 



tail of a comet whose head passed before the beginning of the 
observation. 

We can quantify the geometry of these "cometary" clouds, as 
discussed in the coming sections. Note however, that their basic 
scale depends heavily on the clouds' velocities and on the di- 
mension of the X-ray source, which are both unknown, although 
some constraints can be set. 

The sketch in Fig. [8] summarizes the geometrical constraints 
inferred for the two clouds derived below. This is obtained under 
the assumption that the tail is elongated exactly along the direc- 
tion of motion of the cloud. Moreover, based on the analysis of 
the temporal variation of the covering factor we can only recon- 
struct the shape of the cometary tail projected on the plane of the 
sky. Our data cannot probe any elongation of the tail along the 
line of sight. 

In principle the observed variations can be interpreted with 
other, more complex geometries (e.g. curved clouds with density 
gradients along the direction of motion), but the cometary shape 
is the simplest model that can account for both the variations of 
covering factor and of column density, and which can be inter- 
preted in terms of physically plausible scenarios, as is discussed 
later on in sect. [8] 

4.1. The clouds' velocity 

If we assume that the cloud head is a hemisphere (Fig. [8]), then 
the fact that the initial covering factor of each eclipse is large 
(~65%), but not total, along with the fact that the initial flux 
drop is very sharp, indicate that the cometary cloud head and the 
X-ray source have similar sizes. The sharpness of the flux drop 
at the beginning of each eclipse provides constraints on the ratio 
Rx/Vc between the size of the X-ray source (~ size of the cloud 
head), and the cloud velocity (in the plane of the sky). 

Fig. |9] shows the detail of the 2-5 keV light curve around 
the beginning of eclipse "A", with a binning of 500 sec (top) 
and 100 sec (bottom). The gaps are due to the Earth occultation 
during the Suzaku orbit. The sudden flux drop due to the passage 
of the head of the cloud is marginally resolved, just before the 
Earth occultation gap, and it lasts about one ksec, or less. We 
modelled the expected light curve with a simple model where 
the X-ray source has a radius Rx, the cometary cloud head is 
a hemisphere with radius Rhead ~ Rx, moving with a velocity 
Vc transverse to the line of sight, and covering 65% of the X- 
ray source right after the "head" transit. The minimum possible 
radius for the X-ray source is given by the minimum stable orbit 
in the case of a maximally rotating black hole, which, in the case 
of NGC 1365 implies Rx > 0.5 Rs = 3 10" cm (sect.©. This, 
very conservative, lower limit on the size of the X-ray source, 
along with the short timescale of the beginning of the eclipse (< 
Iksec), sets a lower limit on the velocity of cloud "A", Vc(A) > 
3000 km/s. By using these lower limits (i.e. Rx = 3 10" cm 
and Vc = 3000 km/s) we obtain the light curve shown by the red 
solid line in Fig.|9] A larger size of the X-ray source would imply 
a correspondingly higher cloud velocity, with the constraint that 
Rx/Vc(A) ~ 1 ksec. 

For cloud "B" the eclipse starts within an Earth occultation 
gap. In this case we can only set an upper limit to the ratio 
Rx/Vc(B) < 3 ksec, hence we infer a lower limit on the cloud 
velocity of Vc(B) > 1000 km/s, by assuming for Rx the small- 
est stable orbit for a maximally rotating black hole. 

We note that the lower limit on the velocity of cloud "B" 
is consistent with the average velocity for BLR clouds in- 
ferred from the FNV HM of the broad H/? fine, i.e. 1900 km/s 
(ISchulz et al.L Il999l) . in agreement with the scenario where the 
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Fig. 8. Sketch illustrating the geometry of the cometary clouds that eclipse the X-ray source, as inferred by our observations. For 
the sake of clarity the drawing dimensions are not in the correct scale (in particular the horizontal length of the tails should be much 
longer when compared to the size of the cloud heads and the tail opening angles are exaggerated). The variation of the covering 
factor expected by this model is shown with a dotted line in Fig. [2]:. 
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Fig. 9. Detail of the 2-5 ke V light curve around the beginning 
of eclipse "A". The top panel is with bins of 500 sec. The bot- 
tom panel is with bins of 100 sec and the time axis is further 
expanded around the beginning of eclipse. The red solid line in- 
dicates the modelled light curve where the ratio between the size 
of the X-ray source (~radius of the cloud head) and the cloud 
velocity is Rx/Vc = 10** cm (km/s)"' (see text). If Rx - Rs this 
implies Vc - 6000 km/s. 



eclipsing clouds are part of the BLR, as discussed later on in 
Sect.|6] 

For cloud "A" we infer a velocity larger than expected for the 
average BLR cloud, however not an unre alistic one given tha t the 
power law-like profile of the broad lines lNagao et al.l (l2006h im- 
plies a large fraction of clouds with velocity larger than inferred 
from their FWHM. Yet, if Rx » Rs then the clouds' velocity 
is unrealistically large compared with the observed profile of the 
broad lines. 



4.2. The opening angle of the cometary tail 

The evolution of the covering factor as a function of time pro- 
vides constraints on the geometry of the cloud tail. The model 
used to fit the data (Fig. [8]) depends only on two parameters: 
the initial covering factor (~ 65% both for cloud "A" and "B") 
and the quantity Rx/(Vctan0taii), where ^tau is the opening an- 
gle of the cloud tail. In the case of cloud "B" (where the evo- 
lution of the covering factor is better constrained) we obtain 
Rx/(Vctan0taii) ~ 35.8 ksec. The resulting fit to the data is 
shown with a dotted line in Fig.|2j;. This value, along with the up- 
per limit on the ratio Rx/Vc(B) < 3 ksec inferred above, implies 
an upper limit on the opening angle of the tail of 0taii(B) < 4.8°. 
In the case of cloud "A" the evolution of the covering is de- 
termined more poorly, due to the larger errors. In this case we 
can only set an upper limit to the quantity Rx/(Vctan0taii) > 
47.7 ksec, whose corresponding fit to the data is shown with a 
dotted line in Fig.|2j;. However, as discussed above, in this case 
we have a better constraint on the ratio Rx/Vc(A) ~ 1 ksec, 
hence we can obtain a stronger upper limit on the tail opening 
angle, 0teii(A) < 1.2°. 



4.3. The length of the cometary tail 



As shown in Fig. |2j5, the column density Nh.2 of each cloud 
drops rapidly (the dotted lines represent a simple, empirical 
quadratic fit to the observed evolution of Nh,2). The length 
of the tail is set observationally by our capability to disentan- 
gle the cometary cloud absorption from the (constant) global 
foreground absorber Nh,i, which has a value of 7.5 10^^ cm"^ 
(Table [U. Through simulations we have obtained that our spec- 
tra cannot recover a partial absorber with Nh,2 < 2 10^^ cm"^ 
(i.e. a factor of ~10 below the Nh,2 in the head). With this ob- 
servational "definition", the time elapsed between the passage of 
the cloud head and the time when Nh,2 drops below this limit- 
ing value, along with the lower limit on the cloud velocity in- 
ferred above, provides a lower limit on the length of the tail: 
Ltaii(A) > 2 10'^ cm andLtaii(B) > lO'^ cm. 



5. How common are "comets" in AGNs? 

In this section we discuss whether the cometary geometry of the 
clouds observed by us are peculiar to these observations, or to 
NGC 1365, or whether they may be a more general property of 
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the absorbing clouds in other AGNs, not so far recognized so far 
in the data. 

NGC 1365 may itself appear a peculiar case, since it is the 
AGN where the stro ngest and most frequent Nh vari ations have 
been observed (e.g. iRisaliti et all l2005ai l2007l l2009a) . This is 
partly due to the fact that NGC 1365 is extremely bright in the 
X-ray band, which is a requirement to properly monitor Nh vari- 
ations even on short time scales. Some of the other nearby bright 
AGNs have indeed displayed Nh varia tions on short time scales 
dElvis et al.Ll2004tlPuccetti et all 12007 ). 

However, among the other bright sources, which could in 
principle allow the same investigation, NGC 1365 probably has a 
fortunate vie wing angle. I f the obscuring clouds have a flattened 
distribution (lElvisll2000t IRisaliti et all l2002t iMaiohno et aP. 
1200 U iGaskell et al.'. 7007), coplanar with the dusty circumnu- 
clear absorber (Nenkova et al., 2008), then in most type 1 AGNs 
it is rare to have a cloud transiting along the line of sight. On the 
other hand, in the classical type 2 AGNs, with equatorial viewing 
angle, the number of clouds transiting along the line of sight is so 
large that they make, on average, the absorption nearly constant. 
NGC 1365, which is a type 1.8 AGNs (Schulz et al, 1999), prob- 
ably has an intermediate viewing angle, where there are a signif- 
icant number of clouds transiting the line of sight, but not large 
enough to average out the Nh variations. Therefore, in this sce- 
nario, NGC 1365 is a particularly good laboratory to investigate 
clouds through their transits along the line of sight, even though 
intrinsically it may have no peculiar properties. An alternative 
possibility is that the density of absorbing clouds surrounding 
different AGNs vary from object to object (in this scenario the 
difference between type 2 and type 1 AGN would be mostly due 
to the number of clouds rather than inclination), and NGC 1365 
may have the appropriate density of clouds that makes the proba- 
bility of detecting eclipses high, but not large enough to average 
out Nh variations. 

Assessing the detectability of the cometary geometry of 
the clouds in other objects similar to NGC 1365, as well as 
in other observations of NGC 1365 itself, is not simple. Long 
(> 100 ksec), continuous observations, like the one reported 
here, are required to detect the cometary tail in a transiting cloud. 
This kind of observation is performed rarely. In most cases long 
integration times are obtained through several, discontinuous 
short observations (a few 10 ksec each). Even when such data 
are available the analysis must be tuned to find these types of 
changes. This involves the complex simultaneous fit over the var- 
ious time intervals by leaving free both covering factor and Nh 
for the variable absorber. Luck also plays a role: 1) the transit- 
ing clouds must be moving fast to make the covering factor and 
Nh variations easy to identify in terms of a cometary shape; 2) 
the variation of Nh and of covering factor must be pronounced, 
else the two quantities would result in a degenerate fit; 3) if the 
head of the cloud does not eclipse the X-ray source, and only 
the tail transits in front of the X-ray source, then the evolution 
of the light curve and of the covering factor would appear much 
smoother, and would not be easy to ascribe to the cometary tail. 

To our knowledge, the only previous (nearly) continuous, 
moderately long observations of a bright AGN, with variable ob- 
scuration eve nts, are the XMM lon g observations of NGC 1365 
presented in iRisaliti et al] (l2009ah . iRisa iti et al] (l2 009b) and 
iTurner et al] (l2008h . The observation in iRisahtiet al.. (.2009al) 
does not show an "eclipsing" event, but rather the uncovering 
of the X-ray source due to a "hole" in the clouds distribution, 
hence it is not suitable t o trace t he geometry of the clouds. In 
the data of IRisaliti et all (l2009bh the variation of Nh,2 and CF 
are much milder than observed in our data; as a consequence, in 



iRisahti et al] (l2009bl) it was not possible to leave both parame- 
ters free in all time intervals, since the two quantities would re- 
sult degenerate. Therefore, in Risaliti et al. (2009b) only the cov- 
ering factor was left free to vary among the time intervals and de- 
tails on the structure of the eclipsing cloud could not be inferred. 
The inferred evolution of the CF during cloud transit appears 
relatively symmetric, with no obvious evidence for a cometary 
shape. However, if we had also kept Nh,2 constant in our obser- 
vation, then we would have obtained a more symmetric variation 
of the covering factor: the decreasing Nh.2 obtained by us in the 
second part of each eclipse (Fig.|2]i would be interpreted (with a 
worse ;i^^) as a decreasing covering factor in a model where Nh,2 
is kept constant. Therefore, the eclipse observed in the XMM 
observation may also be cometary-like, but the adopted model 
may have prevented its identification. Yet, the light curve of the 
long XMM observation does not show the sharp flux drop ob- 
served at the beginning of our eclipses "A" and "B". This may 
be either due to the fact that the cloud in the XMM observation 
was moving more slowly (hence the eclipse transition was much 
smoother in time), or to the fact that the head of the cometary 
cloud did not pass in front of the X-ray source (hence only the 
smoother occultation produced by its tail was observed), or more 
simply because the cloud head was "fuzzier". 

Similar c onsid erations apply to the eclipse identified by 
iTurner et al.l pOOS), which was not modelled in terms of both 
variable covering factor and variable Nh, and whose light c urve 
does not show a sharp flux drop, as in IRisaliti et alJ (l2009bl) . 

However, sharp drops of the X-ray flux in few ksec (simi- 
lar to those observed at the beginning of eclipses "A" and "B"), 
have been observed in prev ious observations of other AGNs. 
iMcKernan & Yaqoobid 19981) investigated in detail one such drop 
in an ASCA observation of MCG-6-30-15, which they ascribed 
to the occultation of the X-ray source by a "symmetric" cloud. 
However, they could not investigate the spectral evolution of this 
occultation, which is required to infer the presence of a cometary 
tail. By only using the light curve the beginning of our eclipse 
"A" would be interpreted as an occultation by a nearly symmet- 
ric cloud; only the investigation of the spectral evolution reveals 
the cometary shape. 

Finally, we mention that our team has found evidence for 
cometary clouds also in another source (Risaliti et al. 2010, sub- 
mitt.) through the detailed temporal analysis of X-ray data. 

Summarizing, the cometary shape could be a common prop- 
erty of the X-ray obscuring clouds in general. However, this fea- 
ture is difficult to identify if the observations are not adequate or 
the source not bright enough, so the occurrence rate is presently 
unknown. 



6. Cometary clouds as BLR clouds 

Several papers have discussed that rapid changes (days-hours) of 
the absorbing column density in AGNs support a location of the 
associated ab sorbing clouds within the rad i us of the Broad Line 
Regions ( Elvis et aj ll2004t iPuccetti et all l2007t IRisaliti et all 
2005 a, 2007, 20093B). The rapidly variable absorption in our 
Suzaku data can also be ascribed to BLR clouds. In the following 
we provide a more quantitative justification of this scenario. 

We assume that the absorbing clouds are in Keplerian rota- 
tion around the supermassive black hole. In the case of a spheri- 
cal cloud crossing the line of sight (with the X-ray source smaller 
or comparable to the projected size of the cloud), the distance D 
of the absorber from the black hole can be inferred from the 
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crossing time, traced by the absorption variation, through the 
equation (iRisaliti et al.Ll2002h 



D ~ 6 lO^^ il nl Nh^ Rs [cm] 



(1) 



where td is the crossing time in units of ksec, Ud is the cloud 
density, in units of cm"-', and Nh its column density, in units of 
cm"^. In our case clouds are certainly not spherical. However, 
we can approximate the dense "head" of the cometary clouds 
with hemispheres, as discussed above. In this case td is twice 
the transition interval at the beginning of the eclipses, i.e. 1^ ~ 
4 ksec (by taking the average of clouds A and B, Sect. 14.1b . 
For the column density of the heads we can assume the aver- 
age of those observed at the beginning of eclipses A and B, i.e. 
Nh ~ 1 .5 10^^ cm~^. The gas density of the cloud can be roughly 
estimated by assuming Rhead ^ Rs ~ 6 10" cm (as justified in 
Sect. |4]i, hence n^ - nhead ^ NH/2Rhead ~ 10" cm"-'. In this 
case Eq. [T] gives D < 2 lO'^ c m. This value is consistent with 
that inferred by ( Risaliti et all 12009b,) through an independent 
analysis of other eclipsing events. 

Although this is an approximate estimate, given the vari- 
ous assumptions involved, it is within the radius of the BLR 
inferred for NGC 1365. Indeed, by using the Rblr - Lx rela- 
tion given in'Kaspi et al.' fe OOSi) and the intrinsic X-ray luminos- 
ity inferred in (Risaliti et al.l 12009a, adapted to our new slope). 



L2-iokeV ~ 5 lO'*^ erg s"', we infer a radius of the broad line 
region of Rblr ~ 10'^ cm. This result strongly supports the sce- 
nario identifying the cometary clouds found in the Suzaku data 
with BLR clouds. It should be noted that the BLR radius given 
by the relation in Kaspi et al. (2005) refers to the clouds pre- 
dominantly emitting H/S. However, the radius of the BLR actu- 
ally spans ab out an order of magnitu de as inferred by reverber- 
ation studies (iKaspi & Netzed . [19991) . with high ionization lines 
being emitted predominantly in the inner BLR while low ion- 
ization lines predominantly coming from the outer region. The 
distance inferred for the absorbers (D < 2 10'^ cm), relative 
to the BLR radius inferred for the Hfi emitting clouds, suggests 
that the absorbers observed by us are mostly associated with the 
high-ionization clouds (e.g. those emitting CIV). 

We can infer the mass of individual clouds based on the den- 
sity inferred above. In particular the "head" of the clouds has 

a mass of Mhead ^ nhead mn l^rRhead ~ ^ ^O"" ^0 (^ere we 
have assumed that the cloud head is a full sphere, whose lead- 
ing hemisphere is probed by the observed variations of Nh and 
CF, while the trailing hemisphere is observationally confused 
within the cometary tail). For what concerns the tail, given its 
small opening angle we can approximate Rtaii ~ Rhead ~ Rx S 
Rs = 6 10" cm. We can infer the density of the tail by assum- 
ing that the cross section of the tail is axi-symmetric relative to 



the direction of motion, hence n^ 



ail 



NH,tail/Rtail < 10'" Cm- 



(Fig. |2ll. Out to the length that we can detect trough our ob- 
servations (Figs. |2][8]l we infer a mass of the tail of Mtaii > 
ntaii mn nRl^ Ltaii ~ 10""' Mq. 

From the mass of individual clouds we can derive the mass 
of the BLR if the density of clouds is known. The latter can 
be inferred from the frequency of the eclipsing events. Our 
and past observations of NGC 1365 display eclipses in -1/3 
of the monitoring time, implying that about one hundred of 
absorbing clouds must be present in the volume, "V, probed 
by our Une of sight toward the X-ray emitting source in one 
year, within the region hosting the absorbers. We can write 
^ ~ 2Rhead (Vc ■ lyr) D ~ 5 10""^ 



pc- 



sumed Vc - 2000 km s"' and D ~ 2 lO'^ cm as representative 
of the clouds velocity and of their distance to the center (as dis- 
cussed above). The inferred volume density of the absorbers is 



therefore lOO/'V ~ 2 lO'"* pc ■'. If we assume that this density 
is uniform within the radius of the BLR, then we obtain a total 
number of clouds of about A/c ~ 3 10^. 

The total mass of the BLR inferred from the properties of the 
absorbers is M^^l^ ;« Afc (Mhead + Mtaii) «; 4 10"^ Mq. The latter 
value is five orders of magnitude lower t han inferred from pho - 
toionization models. Indeed, according to 'Bald win et al.l (l2003h . 
the mass of the BLR inferred from photoionization models 



scales as M 



phot 
BLR 



2 10-'*'*(vL,) 



Mq (note that we do 



r; 



not use the monocromatic luminosity as in_ Baldwin et all 
2003). Assuming (vLv)[45o^/L2-iokeV ~ 7 (f^ung et al.. 20l3 
Ric hards etal.Ll2006h . for NGC 1365 we obtain MP^^ ~ 0.7 Mq. 
Accounting for such discrepancy is not easy. One possibility 
is that the X-ray source is much larger than assumed by us 
(Rhead ~ Rx is> Rs). The inferred density of absorbers scales 
as ~ 1/Rx, while the mass of individual clouds scales as ~ R^, 
therefore the mass of the BLR inferred from the absorbers scales 
as ~ Rx. However, to obtain M^^^ similar to M'^^^ would re- 
quire ~ Rx = 200 Rs, not easy to account for by standard mod- 
els. Moreover, as discussed in sect. 14.11 Rx 3> Rs would imply 
clouds' velocities unrealistically high relative to those inferred 
from the broad lines profile. 

Alternatively, a significant fraction of the BLR may not be 
accounted for by the rapidly variable absorbers detected here. 
We may be biased towards the detection of "small clouds", sim- 
ply because the "eclipses" method favors the detection of clouds 
with size similar to the X-ray source; larger clouds may well ex- 
ist, but may be more difficult to detect through Nh variations. In 
particular, part of the "constant" absorber (Nh.i = 7.5 10^^ cm"^) 
may be associated with much larger clouds in the BLR (possibly 
located in the outer part), whose Nh variability is only detected 
on much longer timescales (e.g. see Risaliti et al., 2009a b). The 
mass of such large clouds may account for the bulk of the mass 
of the BLR. 

Yet another possibility is that the mass of the BLR clouds 
inferred by photoionization models may be erroneously (over- 
)estimated. Indeed, such models assume that the photoioniza- 
tion is due to a nuclear point-like source (the accretion disk). 
However, as discussed in the next section, part of the photoion- 
izing flux may be produced locally by the shock generated by 
the interaction between the clouds and the diffuse intracloud 
medium. This is expected to significantly relax the requirements 
of classical BLR photoionization models in terms of mass. 

7. Dynamics and fate of the cometary clouds 

The differential velocity of the gas in the tail relative to the head 
(Vhead - Vtaii) is simply given by the sound speed of the gas in the 
cloud (~ 10 km/s), since behind the head there is no hot intra- 
cloud medium ("cleaned" by the passage of the head), nor other 
agents preventing the gas from expanding freely. Essentially, the 
tail follows the head in its motion, and lags behind only by a 
small fraction of the bulk velocity. 

We can estimate the cloud mass-loss through the tail, which 
is given by 



Mhead = (Vhead " Vtail) ntail mn TtR. 



tail 



(2) 



We obtain a mass-loss rate of about Mhead ~ 3 10" Mq yr" . By 
comparison with the mass of the clouds head inferred in the pre- 
vious section we infer a lifetime of the clouds of thfe ~ 2 months, 
i.e. shorter than their orbital period (lyr). A similar conclu- 
sion on the short hfetime of BLR clouds was reached also by 
iMathews & Ferlandl (Il987h based on the expected interaction of 
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BLR clouds with the warm intracloud medium (this issue will be 
discussed further in the next section). 

The latter constraint has important implications. The BLR 
in most type 1 AGNs is observed to be present (and often nearly 
constant) for more than 10-20 years. As a consequence, the much 
shorter lifetime of the BLR clouds implies that they have to be 
continuously replenished, and at a very small radius (< 10'^ cm). 
Direct infall at such small radii is unlikely (inflow must occur 
through the accretion disk). The more likely possibly is that BLR 
clouds are continuously produced by the accretion disk, as sug- 
gested by the model in Elvis (2000). 

It should be noted that tufe depends on our assumption on Rx, 
since Mhead/Mhead depends linearly on Rhead ~ Rx- However, 
to make the lifetime of the absorbing clouds longer than sev- 
eral years (to account for the observationally continuous pres- 
ence and "stability" of the BLR) would require Rx » Rs, 
which runs into other physical problems discussed in sect. |6l 
Yet, within this context, we note that the argument on the short 
timescale only applies to the "small" cometary clouds detected 
through the short timescale Nh variability. The putative larger 
BLR clouds mentioned in the previous section may not suffer 
the same timescale problem if they are much more massive. 

The global mass loss of BLR comets is Mtot - A/cMc ~ 
10"^ Mq yr~'. Modelling the fate of such gas lost by the comets 
is very difficult and goes beyond the scope of this paper. Here we 
only mention two possibilities. If the density drops significantly, 
the eff'ective ionization parameter may become high enough to 
make the gas thermally unstable and heating up to 10^ K, hence 
contributing to the hot intracloud medium (see next section). 
The gas lost by the cometary clouds may also be accelerated 
by the radiation pressure contributing to a fraction of the wind 
often observed in Seyfert nuclei, and who se outflow rate is gen- 
erally estimat ed to be around 0. 1 Mq vr~' (iKrongold et alll2007l: 
ICrenshaw et a l.. 2003; Andrade- Vel azquez et 
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8. The nature of the cometary clouds 

It is beyond the scope of the present work to define a detailed 
physical model accounting for the nature of the "cometary-like" 
structure of the clouds. In this section we only try to identify the 
most plausible scenarios, by exploiting the constraints provided 
by the observations discussed in the previous sections. As dis- 
cussed above, these cometary clouds are most likely part of the 
BLR. Therefore, any discussion on their nature is necessarily in- 
terlaced with the nature o f the BLR clou ds in general, which is 
still matter of debate (see iNetzeii 120081 for a review), although 
here we mostly focus on their cometary shape. 

First, although we use the term "comets", these are not 
comets in the solar system sense, i.e. their tail is not due to radi- 
ation pressure on the gas. Indeed, the latter would make the tail 
elongated radially and yielding symmetric variations of the Nh 
and of the covering factor, in contrast to what is observed. Yet, 
as mentioned in sect. |4] we cannot exclude a component of the 
tail elongated along our line of sight. 

In Sect. |6] we inferred that each cloud head must have a 
gas density of about lO" cm"^, while the density drops by 
about an order of magnitude along the tail. The large den- 
sity of the head may suggest that we are actually seeing 
the (bloated) atmosphere of a star passing along the line of 
sight. Models ascribing the BLR clouds to the atmosphere 
of ( bloated or evolved) stars have been proposed in the past 
JScoville & Norman, 1988 ; Alexander & N etzer, 1994, 199^ 
and it was also suggested that such s t ellar a tmospheres may pro- 
duce contrails dScoville & Normanl Il995[) . which may mimic 



the "tails" observed in our clouds. However, as discussed in 
sect.|6l the inferred density of the absorbing clouds (~ lO''* pc"^) 
is by far larger than the stellar density in the center of galaxies. In 
the light of the present X-ray observations, modelling the BLR 
structure in terms of stellar atmospheres seems not to be a phys- 
ically viable interpretation. The absorbing (and BLR) clouds are 
likely purely gaseous (non self-gravitating) objects. Yet, we can- 
not exclude that some of the eclipsing events are due to transiting 
stellar atmospheres. 

Within the general framework of AGN central regions, these 
"cometary-clouds" should be moving through a diffuse, (rela- 
tive ly) hot in tracloud medium (HIM), as originally discussed in 
Krohk et alJ (11981) . This scenario has been subject to revisions 
and objections and, at present, the tempe rature of such diffuse 
medium is estimated to be Thim ~ 10^ K (iPerrv & Dv son'.' 19851 



K (iPerrv & 



Netzer, 199^ 



Fabian etal, 1986; Collin-Souffrinet 
Kroli^ Il999i) . This intracloud medium might well be in out 
flow, but it is highly unlikely that such medium is in orbital 
motion around the central black hole, so that we can exclude 
that the absorbing clouds with the inferred velocities are co- 
moving with this ambient medium (at least for what concerns 
the orbital component of their motion). A simple estimate of the 
HIM sound speed shows that these cometary clouds must indeed 
be moving supersonically through the HIM. In fact, the sound 
speed in the HIM is (cJhim ~ 290 (Thim/ 10^)'^^ km s"', while 
the inferred velocities of the "cometary" clouds are significantly 
larger (Sect. l4Tb . The inferred Mach numbers are Ma > 20 and 
Mb > 7 (recall that only lower limits on the clouds velocities 
were obtained in Sect. 14. lb . These imply Mach cones with open- 
ing angles (< 1.5° and < 4° for clouds A and B, respectively) 
consistent with the tails opening angles inferred by us through 
the variation of the covering factor 

The supersonic motion of dense clouds in a less dense 
medium has been extensively studied in the pas t, although 
mostly in very different astrophysical contexts (iFalle et al.L 
l2002l:|Pittard et al.Ll2005alfbh . Such theoretical studies show that 
in general a sort of ("cometary") tail, behind the front shock does 
indeed tend to form. In numerical simulations the specific re- 
sult (shape, physical properties, etc..) depends on the particular 
physical processes accounted for (for example, whether radia- 
tive cooling is included or not) and on the physical parameters, 
especially on the density contrast between unshocked ambient 
medium and cloud material, and on the Mach number. In any 
case, a system of front shocks (bow-shock like) at the head of 
the cloud, and a sort of Mach cone, prolonging the bow-shock at 
the rear of the cloud-head itself, must form. The "cometary tail" 
is made up by the debris of trailing material lost by the cloud 
head, because of the action both of Rayleigh-Taylor and Kelvin- 
Helmholtz instabilities (Kroliket al., 1981; Mathews & Ferlandl 



1987) and of oth e r possi ble disrupting mass loss mechanisms 



(e.g. iPittard et al.L 120051 ^, and reference therein). Moreover, as 
discussed above, the lack of a confining medium on the trailing 
side of the cloud head naturally leads to gas loss at the sound 
speed (~ lOkm/s) due to free expansion, even in absence of any 
hydrodynamical instability. An elongated structure, trailing be- 
hind the cloud denser head, is therefore expected by models, and 
also with a cometary-shape. However, a more quantitative com- 
parison with our data is extremely complex, since it requires that 
such models are modified to match the physical conditions ob- 
served in AGNs and in NGC 1365 in particular 

It is important to note that, as first pointed out by 
ICollin-Sou ffrin et al. ( 1988), the shock generated by the super- 
sonic motion of the BLR clouds into the HIM is expected to gen- 
erate UV and X-ray radiation that can contribute significantly to 
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the ionization in situ of the cloud itself. This effect requires a 
substantial revision of the classical photoionization models of 
the BLR, which assume photoionization by a nuclear, point-like 
radiation source. 

Finally, we mention that magnetic fields can provide an addi- 
tional mechanism that could elongate the circumnuclear clouds 
into a filamentary, and possibly cometary, shape. Strong mag- 
netic fields are expected to surround accreting black holes, and 
have been pr oposed as an alternative mechanism to confine the 
BLR clouds ( Reei [19871: fPerland & ReeslfToSSl) . with the effect 
of making them elongated along the field. Inferring the actual 
shape of the clouds expected to originate from such a mecha- 
nism is also very complex and beyond the scope of this paper, 
a more detailed analysis of these effects is deferred to a future 
paper. Here we only note that the past models of magnetically 
confined BLR clouds do predict elongated structures, but do not 
obviously predict dense heads as those observed here. 



9. Summary 

We have presented a detailed spectral analysis of a long (~ 
300 ksec), continuous Suzaku observation of the Seyfert nucleus 
in the galaxy NGC 1365. The spectrum shows evidence for an 
absorption component that is variable in time, both in terms of 
column density and in terms of covering factor, which is ascribed 
to clouds eclipsing the X-ray source. This is the first time that a 
temporally resolved X-ray spectral analysis is able to break the 
degeneracy between the evolution of the column density (Nh) 
and the covering factor (CF) of the X-ray absorber. 

We identify two main eclipses. The temporal evolution of 
each eclipse is far from being symmetrical. The initial occulta- 
tion is very rapid (within ~ 1 ksec) and covering about 65% of 
the source. Subsequently the covering factor increases, but more 
slowly, reaching unity in about 50 ksec. The absorbing column 
density of the cloud is about lO^"* cm"^ at the beginning of each 
eclipse and decreases afterwards. These results are inconsistent 
with a spherical geometry for the absorbing clouds. The most 
likely geometry compatible with the observations is an elon- 
gated, "cometary" shape, with a dense "head" (n ~ 10" cm"-', 
consistent with that expected for the clouds of the BLR) and with 
a dissolving and expanding tail. 

The data allow us to quantitatively constrain the geometry, 
dynamics and location of such cometary clouds. The cometary 
clouds are probably located at a distance of about 2 10'^ cm from 
the nuclear black hole, well within the estimated BLR radius (~ 
10'^ cm), strongly supporting the association of these absorbing 
systems with BLR clouds (at least with the inner, high-ionization 
ones). The cometary clouds "head" must have a size comparable 
to the X-ray source and must be moving with velocity higher 
than about 1000 km/s (consistent with the velocity expected for 
the BLR clouds). The cometary tail must be longer than a few 
times 10'^ cm. The tail opening angle must be very narrow, less 
than a few degrees, and consistent with the opening angle of the 
Mach cone expected from the supersonic motion of the cloud 
into the hot intracloud medium. 

We suggest that such cometary clouds may be common to 
most AGNs, but have been difficult to recognize in previous X- 
ray observations. 

We estimate that the cloud "head" loses a significant frac- 
tion of its mass through the cometary tail, which is expected to 
cause the total cloud destruction within a few months. If these 
clouds are representative of most BLR clouds (or at least the 
high-ionization ones), our result implies that the BLR region 



must be continuously replenished with gas clouds, possibly from 
the accretion disk. 

We briefly discussed the possible nature of such cometary 
clouds. The most likely scenario is that the tail is made of 
gas lost by the cloud head through hydrodynamical instabilities 
generated by its supersonic motion through the hot intracloud 
medium. 

We can estimate the mass of the absorbing clouds (Mdoud ~ 
10""^ Mq) and their total number within the central region 
(Nc ~ 3 10^). The inferred total mass of the BLR is about 
4 10"^ Mq, which is two orders of magnitude lower than the 
BLR mass inferred from photoionization models. The discrep- 
ancy may originate from a population of large, massive BLR 
clouds not identified in our eclipsing studies. Alternatively, pho- 
toionization models may overestimate the BLR mass. In particu- 
lar, UV and X-ray radiation produced by the shocks generated by 
the supersonic motion of the clouds may provide a local source 
of ionizing photons, not accounted for by classical photoioniza- 
tion models that assume a central point-like radiation source. 
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